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ABSTRACT 


The Chandra High Energy Transmission Gratings (HETG) Orion Legacy Project (HOLP) is the first compre¬ 
hensive set of observations of a very young massive stellar cluster which provides high resolution X-ray spectra 
of very young stars over a wide mass range (0.7 - 2.3 Mq ). In this paper, we focus on the six brightest X-ray 
sources with T Tauri stellar counterparts which are well-characterized at optical and infra-red wavelengths. All 
stars show column densities which are substantially smaller than expected from optical extinction indicating that 
the sources are located on the near side of the cluster with respect to the observer as well as that these stars are 
embedded in more dusty environments. Stellar X-ray luminosities are well above lO*^* erg s“* , in some cases 
exceeding lO*^^ erg s“* for a substantial amount of time. The stars during these observations show no flares but 
are persistently bright. The spectra can be well fit with two temperature plasma components of 10 MK and 40 
MK, of which the latter dominates the flux by a ratio 6:1 on average. The total EMs range between 3 - 8x lO^'* 
cm“^ and are comparable to active coronal sources. Limits on the forbidden to inter-combination line ratios in the 
He-Like K-shell lines show that we observe a predominantely optically thin plasma with electron densities below 
10*^ cm“^. Observed abundances compare well with active coronal sources underlying the coronal nature of these 
sources. The surface flux in this sample of 0.6 to 2.3 Mq classical T Tauri stars shows that coronal activity and 
possibly coronal loop size increase significantly between ages 0.1 to 10 Myrs. 


1. INTRODUCTION 

The Orion Nebula Cluster (ONC) is a complex stellar for¬ 
mation region that hosts a variety of young stellar objects 
in terms of mass, age, configuration and evolutionary stages. 
The cluster is part of the Orion A molecular cloud which 
is hos t to a hierarchical structure of ongo ing star formation 
cells (iBallv, O’Dell & McCaughreanll2000h . The part of the 
ONC we generally refer to is a somewhat older formation bub¬ 
ble located at the foreground of the main molecular cloud. Its 
massive stars as members of the Orion Trapezium - 0* Ori 
C and 9^ Ori A - are the main sources of illumination and 
ionization of the Orion Nebula (M42). The ONC hosts one 
of the largest assembly of young stars in Orion with about 
80% of its members being younger than 1 Myrs. There are 
over 3000 stars in the immediate vicinity of the Orion Trapez¬ 
ium leading to an average stellar density of about 250 stars 
per pc^ within a radius of about 3 pc. About 1600 stars are 
optically observed and have been classified to some limited 
extent through spe ctroscopic and photometric measurements 
dHillenbrandll 1 997h . O ver 2000 stars have be en observed in the 
IR band with 2MASS (ISkrutskie et al.ll2006h and gr ound based 
surveys (iMuench et al.ll2002HRobberto et al.ll2010h . 

X-rays from Orion a nd specifically the O NC were first 
discovered with Uhuru (iGiacconi et al.l Il972h as the bright 
X-ray source 3U0527-05. Half a decade later observa¬ 
tions with the ANS satellite suggested a more extended 
emission region and it was first suggested that the X- 
rays are emission from coronae around T Tauri stars 
(Iden Boggende et fflll978h. This suggestion was finally con¬ 
firmed b y the Einstein dEeigelson & Decamplil Il98lh and 
ROSAT dGagne, Caillault & StauffeJll995h X-rav Observato¬ 
ries. But it was the Chandra X-ray Observatory launched in 
1999 that provided the bulk of our current knowledge of X- 


ray detections, identifications, and basic spe ctral properties. 
The C handra Orion Ultradeep Project (COUP, iFeigelson et al.l 
l2005h detected 1600 X-ray stars and measured column densi¬ 
ties, source fluxes, and basic X -ray spectral and photometric 
parameters dGetman et al.ll2005h . 

X-ray emissions from young pre-main sequence (PMS) stars 
with masses 0.1 Mq to about 2 Mq have 10^ times the flux than 
stars on the main sequence. The bulk of emissions in PMS stars 
is in the optical band. The ratio of X-ray to bolometric lumi¬ 
nosity in these stars lies between 10“"* and 10“^, close or at the 
saturation threshold. Studies of the X-ray emissions are specif¬ 
ically sensitive to coronal activity and to some extent - in con¬ 
junction with UV emissions - to stellar surface shocks from ac¬ 
cretion. In fact, studies within the last two decades revealed that 
while these X-rays are primarily attributed to enhanced coro¬ 
nal activity of the star itself, many of these young PMS stars 
also generate X-rays through accretion from a p r otostellar disk 
(Kastner_et^ 200^ 20041 Schr nitt et al.l l2005l iGiinther et aH 
l2006t Argiroffi, Maggio & PeresE007 ). The cases of TW Hya, 
BP Tau, V4046 Sgr, and MP Mus - all are known to be older 
classical T Tauri stars (CTTS) - are indeed quite enigmatic be¬ 
cause here the bulk of the X-ray emission is observed below 
1.5 keV (above 8 A) with peculiar line ratios indicative of ac¬ 
cretion shock emissions. CTTS by definition exhibit strong hy¬ 
drogen Balmer (Ha) emissions and show stong IR excess emis¬ 
sions up to 100 microns indicating the presence of active ac¬ 
cretion disks. X-ray studies of the Taurus regi on also indicate 
that accretion signatures are present in CTTS (iTelleschi et al.l 
l2007l lAudard et al.ll200'^ . The stars mentioned above are the 
best studied CTTS in X-rays today. In contrast, most PMS 
stars, CTTS and older weak-lined T Tauri stars (WTTS), are X- 


bright due to active coronae. iKastner et al 


( WTTS) . 
i1 (I2004II 


clearly 


ray 

showed that coronal emissions are primar ily responsible for 
the WTTS HID98890, |Telleschi et al.((l200^ showed that many 
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CTTS have hard spectra with substantial emissions up to 10 
keV, which are definitely not due to accretion shocks. 

Evidence that most CTTS have X-ray spectra dominated by 
hard emission and produced by coronal activity comes from 
very young stellar clusters and specifically the ONC as the 
closest massive cluster to the Sun. Chandra observations of 
many young clusters in the stellar neighborhood revealed thou¬ 
sands of X-ray sources, mos t of them CTTS, with X-ray spec- 
tra as hard as 2 - 3 keV ( FeigeEonet^ 2005L IWolk et ^ 
2001 iRho etaklllOOl iTownslev et alJ 1201 Ih . In the ONC 


Preibisch et alJ (l2005h showed that in a COUP sample of 600 


X-ray sources, which are reliably identified with optically well 
characterized T Tauri stars, the plasma temperatures obtained 
from the X-ray spectra not only are much hotter than usually 
observed in main sequence stars but also show correlations re¬ 
lated to coronal activity. Most details of coronal properties in 
PMS stars are still unknown as well as its evolution towards a 
much quieter state when these stars reach the main sequence. 

In this paper we present an in-depth analysis of the X-ray 
properties of six classical T Tauri stars in the ONC for which 
we have sufficient spectral data. These stars are some of the 
bri ghtest and most active CTTS in the COUP sample described 
bv fPreibisch et alJ (l2005h . 


2. OBSERVATIONS AND DATA REDUCTION 

The Orion Trapezium was observed with the Chan¬ 
dra hi gh energy transmi ssion grating spectrocmeter 
HETGS (ICanizares et al.ll2005h seyenteen times between 1999 
and 2008. Table [T] summarizes the obserying parameters and 
exposures. While most obseryations were performed by point¬ 
ing towards the brightest Orion star 0' Ori C (RA: 05:35:16.46, 
Dec: -05:23:22.85), some pointings were up to 1.5’ off that di¬ 
rection. Since the six PMS stars we focus on here are within 40" 
of 0' Ori C, this created some yery unfayorable combinations 
of off-axis angles and roll angle oyerlaps. Specifically OBSIDs 
4473 and 4474 are seyerely affected and we omit the latter from 
the analysis. In OBSID 4473 we only haye useful spectra for 
VI 399 Ori. Note, the count rate in Table [T] only shows the 
cumulatiye rates of all six sources after spectral cleaning (see 
below). While the total exposure of all obseryations is about 
585 ks, at the end only 250 ks on ayerage are ayailable for each 
of the six stars in the sample analyzed here. 

All obseryations were reprocessed using CIA04.6 with the 
most recent Cl AO CALDB products. We used standard waye- 
length redistribution matrix files (RME) and generated effectiye 
areas (AREs) using the proyided aspect solutions 0. Note that 
for the HETGS spectra the RME is fairly independent of focal 
plane temperature, howeyer the order sorting has to be adjusted 
to the different pulse height distributions of the CCD at -110 
C in OBSIDs 3 and 4. Here we followed the steps applied in 
the data reduction for the Trape zium stars described in pre yious 
analyses using these data sets (ISchulz et al.ll2000l l2006l) . Eor 
all the HETGS obseryations we generated spectra and analysis 
products for the medium energy gratings (MEG) H-1 and -1 or¬ 
ders, as well as for the high energy gratings (HEG) H-1 and -1 
orders. 

^ see http: //asc. harvard. edu/ciao/threads / 



Fig. 1.— Plots of grating dispersion (tg_mlam) versus CCD array pulse 
height channels (energy) for one OBSID for MEG (left) and HEG (right). The 
source dispersion of the bright star MT Ori (red) overlaps significantly with 
another source (black) rendering the MEG portion unuseable. 

A crowded field of stars and obseryations at yarious differ¬ 
ent roll and off-axis angles create additional challenges in the 
proper spectral data extraction. In order to clean the data there 
are seyeral steps to conside r. These are described in detail in 
[Huenemoerder et aP (l2009l) for the extraction of the young and 
bright intermediate Orion Trapezium binary 0' Ori E. The first 
one concerns the yarious off-axis angles of the sources with 
respect to the telescopes optical axis. The HETG spectral res¬ 
olution degrades significantly once a source is well beyond 1 ’ 
from the optical axis of the telescope. 

As a next step the orders from the grating dispersion are 
sorted by the pulse heights (PHA) of the CCDs. Howeyer 
sources with rolls causing oyerlaps of their dispersed spectra 
from other sources spatially will suffer from sharing PHAs of 
dispersed orders from these sources. The result cannot be re- 
solyed and these data haye to be discarded. Depending on roll 
angle, single grating arms are affected differently and each arm 
for each source has to be inyestigated separately. Note that these 
oyerlaps can be caused by any of the sources in the field includ¬ 
ing the brig ht sources of the Orion Trapezium (ISchulz et al.l 
1200 iL 1200^ . Pigure[T] shows an example of such a case in MT 
Ori data where PHA orders of two different sources oyerlap. In 
all of these cases the grating orders cannot be used. On ayerage 
these occurrences remoye about 40% of the ayailable exposure 
for each source. 

In a separate step coincidences of other source zero orders 
located on HETG dispersion arms haye to be eliminated as they 
appear as spurious line emissions in the dispersion. Many of 
these coincidences can be allocated using the source locations 
from iGetman et al] (l2005h . Howeyer, since there are consider¬ 
able numbers of X-ray flaring sources in the ONC, a final yisual 
inspection of each dispersion arm was warranted. In this respect 
any line feature in the four grating arms for each source had to 
be consistently present, otherwise we suppressed that particular 
bandpass in the dispersion. The rates in Table [1] are calculated 
from the dispersed spectra after all cleaning steps were com¬ 
pleted. 


3. GENERAL SOURCE PROPERTIES 


3.1. Optical Properties 

Physical parameters of yery young PMS stars are hard to de¬ 
termine as at that young stage much of the central star is still 
enshrouded in natal circumdisk material. This means that ac¬ 
curate determinations of extinction and accretion luminosities 
are difficult. A summary of possible biases in the age determi¬ 
nation of young PMS stars can be found in ISoderblom et al.l 
(I 2 OI 4 I 1 . The standard for ONC members was set by a study by 
iHillenbrai^ ( 1997h. which wa s exclusiyely used in the COUP 
analysis (iPreibisch et al.l2()0^ . A more recent improyed multi- 
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Table 1 

Observation Log 


ObsID 

Start Date 
[UT] 

Start Time 
[h:m;s] 

Exposure 

[ks] 

rate^ 
10 “' c/s 

3 

1999 Oct 31 

05:46:18 

49.482 

3.9 

4 

1999 Nov 24 

05:36:51 

30.914 

1.5 

2567 

2001 Dec 28 

12:24:53 

46.357 

2.9 

2568 

2002 Eeb 19 

20:28:38 

46.334 

2.4 

4473 

2004 Nov 3 

01:47:00 

49.119 

0.3 

7407 

2006 Dec 3 

19:06:43 

24.633 

1.6 

7410 

2006 Dec 6 

12:10:32 

13.067 

2.4 

7408 

2006 Dec 19 

14:16:25 

24.862 

2.4 

7409 

2006 Dec 23 

00:46:36 

27.086 

2.1 

7411 

2007 Jul 27 

20:40:17 

24.631 

0.7 

7412 

2007 Jul 28 

06:15:04 

24.834 

0.9 

8568 

2007 Aug 6 

06:53:03 

35.860 

1.1 

8589 

2007 Aug 8 

21:29:30 

50.402 

1.5 

8897 

2007 Nov 15 

10 :02:11 

23.644 

1.3 

8896 

2007 Nov 30 

21:57:29 

22.657 

1.3 

8895 

2007 Dec 7 

03:13:02 

24.851 

1.5 


f sum of 1st order dispersed count rates for all six sources. 
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Fig. 2.— These lightcurves were compiled from HETG 1st order fluxes where available. We abstained from using zero order data because of pile up issues. 
Observations are separated by the dotted lines. The UT times are listed in Tab.[T] 


color optical survey provided updated extinction values and ef¬ 
fectiv e temperat ures for a variety of ONC sta rs dPa Rio et al.l 
2010l) (see also [Herczeg & Hillenbrandll2014l) . The study by 
Da Rio et alJ j2010l) produced a new Hertzsprung-Russell dia- 
gram of the ONC popula t ion. U sing the tracks calculated by 
Siess. Dufour & Eorestinil (|2000|) . this study provided new val¬ 
ues for physical parameters such as Ai/, accretion luminosity 
fractions 'LacJ^tot, masses M* and stellar age ?istar- The values 
for the six stars in our sample are summarized in Table |2] Eor 


the stellar ages we also include the iHillenbrandl (Il997h results, 
a;, 97 , for reference purposes. 

The effective temperatures for all six stars range from ^ 4000 
K (LQ Ori) to about 5300 K (V348 Ori), which puts these stars 
currently into the G/K type category. Except for LQ Ori, the 
stars still produce a significant amount of luminosity from ac¬ 
tive accretion with fractions between ^ 20 and 40 %. Masses 
and ages have ben determined from PMS tracks provided by 
the model from ISiess. Dufour & Eorestinil (l2000t) . which gives 
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Table 2 

Low-Mass PMS Stars: Optical Properties 


Name 

log Teff 

[K] 

log U 
Lq 

AV 

log Lflcf/L* 

M* 

Mq 

R* 

Rq 

log a* 
[yr] 

log a ;,97 

[yr] 

MTOri 

3.66 

1.08 

1.89 

-0.76 

1.99 

8.2 

5.6 

6.2 

LQOri 

3.60 

0.82 

0.46 

-2.42 

0.70 

7.7 

4.8 

6.0 

Pari 842 

3.75 

0.64 

2.84 

-0.37 

1.56 

3.3 

7.0 

6.4 

V13990ri 

3.71 

0.81 

1.61 

-0.74 

2.28 

4.6 

6.4 

6.4 

VI229 Ori 

3.72 

0.86 

2.98 

-0.45 

2.22 

5.0 

6.5 

6.5 

V348 Ori 

3.72 

0.91 

2.97 

- 0.68 

2.33 

4.7 

6.5 

6.3 


a mass range between 0.7 Mq and 2.3 Mq and ages between 
7x 10"^ yr and 1x10^ yr. Most accurate knowledge of these two 
quantities is essential because they provide the stellar radii we 
need to interprete the coronal emissions results. 

3.2. X-Ray Properties from the COUP Study 

The sources in this sample are the X-ray brightest stars within 
a 1.5’ radius around the Orion Trapezium. The source selection 
is purely based on their presence in the optimal field of view 
(FOV) of the Chandra HETG spectrometer. The X-ray prop¬ 
erties of these stars as determined from the COUP dataset are 
lis ted in Table and m ore o f their global propertk s are shown 
in iGetman et aLl ( l2005l) and iPreibisch et ai.l (l2005l) . The COUP 
properties of our sample are affected by photon pileup up to 
100 % and in all the cases mitigating steps had to be taken in the 
COUP analysis. Absorption columns were found to be signifi¬ 
cantly lower than in the bulk of X-ray sources ranging between 
1X10^° cm“^ and 2x 10^* cm“^. The spectra were fitted by a 
two temperture plasma model yielding average temperatures of 
8 x 10^ K and 1.9 X 10^ K, up to over a factor t wo lower than 
most other COUP souces (IPreibisch et al.ll200^ . Emissivities 
are 2.6x 10^"^ cm“^ and 5.5 x 10^"* cm“^ and comparable to the 
entire COUP sample. 

3.3. X-Ray Fluxes from HETG 

We fit the spectra with a two component plasma model 
using the Astroph ysical Plasma Emission Database APED Q 
dFoster et al.ll2012l) . In this first fit we do not constrain parame¬ 
ters such as abundance or absorption, but model the continuum 
as accurately as possible. This first fit already provides reliable 
spectral fluxes and luminosities. The values shown in Table |4] 
are the stellar surface fluxes and surface integrated luminosities 
for a distance of 450 pc. We use the latter distance in order to 
directly compare the results to published COUP data. The ac¬ 
tual dis tance toward the ONC has re cently been determine d to 
412 pc (iReid et al.l2009t) . As done in IPreibisch et al.l(l2005h the 
stellar X-ray flux is defined as a surface flux and is the luminos¬ 
ity divided by the stellar surface area. 

The most luminous PMS star in the sample is MT Ori with 
8.4x 10^' erg s“' , the least luminous star is 3.1 x 10^' erg s“' , 
but all luminosities are lower than deduced in the COUP anal¬ 
ysis. The ratio of X-ray to bolometric luminosity is log Lx/Lhoi 
= -2.77±0.08 and very similar for all sources. Detected X- 
ray fluxes range between 1x10“'^ erg cm“^ s“' and 3x10“'^ 
erg cm“^ s“' . For the luminosities we corrected these fluxes 
for the absorption column. The stellar X-ray surface fluxes (F^: 
in Tab. HI were then determined by dividing these luminosi¬ 


ties by the stellar surface area using s tellar radii associated with 
the mo del parameters in Tablel^and lSiess, Dufour & Forestinil 
(I 2 OOOI) tracks. 


3.4. Variability 

In order to determine source variability we use only the rate 
deduced from the first order HETG spectra in order to avoid 
any biases from pileup in the zero order point sources. The 
light curves are shown in Fig.|2l The 1st order count rates are 
plotted for all observation segments, which in the light curves 
are separated by dotted lines. The observing times of each seg¬ 
ment thus appear on a continuous scale even though the seg¬ 
ments are weeks, months, or even years apart. Some segments 
do not show data, this does not mean that there is no source 
flux but merely that we do not have valid data coverage after 
the confusion analysis. The total observing time is 5.85 x 10^ s. 

We do not observe any significant flaring activity in each of 
the observing segments. The sources appear persistent in most 
observing segments with limited short term flux variations. No¬ 
table rate variations of a few factors on time scales of a few 
hours appear in LQ Ori and VI399 Ori. The larger changes in 
flux appear in between observations. The largest change be¬ 
tween observations appears in VI399 Ori with almost a factor 
10 difference in mean flux per observation. There might be flare 
onsets in the the first segments of Par 1842 and V348 Ori, but 
these appear at the end of the observations and are thus hard to 
classify as a flare. Note that even in the COUP sampl e these 
sources do not have a flare history (iGetman et al.ll200^ . 


4. SPECTRAL ANALYSIS 

The modeling of the spectra and the X-ray line emission is 
done in several steps. The models used in the analysis are two 
temperature APED model spectra for the overall spectral anal¬ 
ysis and local continua and Gaussian functions for the emission 
line analysis. A two temperature model became necessary and 
sufficient since a one temperature unconstrained fit produced 
strong residuals and some unfeasible abundance parameters and 
models with more temperature components did not appear sta¬ 
tistically significant. Hence the two temperature APED models 
are the most consistent way to fit the spectra of all six sources. 
Spectra and fit models are shown in Eigs.[3]and |4] 


'* http://www.atomdb.org ISmith et alJ2001h 
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Fig. 3.— Two temperature APED fits of the combined 1st order HETG spectra of for the brightest source MT Ori and the weakest source Par 1842. 
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Fig. 4.— Unfolded two temperature APED fits of the combined 1st order 
FIETG spectra of four of the Orion sources. 


4.1. Abundance Line Ratios 

Most of the power of the HETG spectra comes from the de¬ 
tection of distinct X-ray line emission. The measurement of K- 
shell line properties puts significant constraints on the physical 
nature of the emitting plasma. Collisional plasmas at temper¬ 
atures above 10® K emit K-shell lines from all major cosmic 
abundant elements in the available wavelength range from 1.7 


to 25 A of the HETG spectrometer. In each spectrum we de¬ 
tect over 100 X-ray lines of which about 70% are very weak Fe 
lines. In all the spectra the lines a ppear unresolved and j match 
the expected locations (see also [Huenemoerder et aP (l2003h 
with no apparent shifts and line broadenings. 

The HETG band covers most abundant non-iron H-like lines 
and He-like line triplets except for C and in most cases N. Mea¬ 
sured line fluxes are shown in Table|5] Many lines, specifically 
for the He-like triplets appear spurious or we only see upper 
limits. The upper limit line flux in our sample is about 2x 10“’ 
photons cm“’ s“'. For high Z He-like triplets such as Ar, S, and 
Si we only see a blend and significantly model the resonance 
and forbidden line components. The detection of higher Z lines 
from Ca and Fe are more rare. 

Ratios of the line fluxes are used to determine a common 
abundance distribution with a first goal to constrain the APED 
plasma model fits. Here we utilize temperature insensitiv e 
abundance ratios as u s ed by i Drake. Testa & Hartmanril (l2005h : 
iDrake & Testa J 2005 Liefke et al.ll2008l) . and further devel- 
oped bv iHuenemoerder et ^ ( 20131) . which read as 


Ai Fiff + aiFif]er ,,, 

A j Fjfl + a2Fjfler 

where are abundances of element i and j, FiH,FjH 

the line fluxes (in photons cm“’ s“') of the correspond¬ 
ing H-like ions, FiHer,Fj^Her the resonance line components 
of the corresponding He-like line triplets. The coeffl- 


Table 3 

Low-Mass PMS Stars: COUP Properties 


Name 

COUP# 

Spec. Type 

RA 

Dec 

logN/, 

(cm-’) 

kT, 

(keV) 

kTz 

(keV) 

log EMi 
(cm-’) 

log EM 2 
(cm-’) 

log 

(erg s-‘ ) 

MT Ori 

932 

K2-K4 

05 35 17.94 

-05 22 45.5 

21.17 

0.76 

1.63 

54.8 

55.1 

32.18 

LQ Ori 

394 

K2V 

05 35 10.73 

-05 23 44.6 

20.00 

0.68 

1.71 

54.0 

54.1 

31.32 

Par 1842 

689 

G7-G8 

05 35 15.23 

-05 22 55.7 

20.85 

0.64 

1.64 

54.2 

54.8 

31.82 

VI399 Ori 

1130 

G8-K0 

05 35 21.04 

-05 23 49.0 

21.28 

0.71 

1.51 

54.4 

54.4 

31.67 

VI229 Ori 

965 

G8-K0 

05 35 18.35 

-05 22 37.4 

21.25 

0.64 

1.63 

54.4 

54.8 

31.89 

V348 Ori 

724 

G8-K0 

05 35 15.62 

-05 22 56.4 

20.46 

0.80 

1.45 

54.3 

54.6 

31.73 
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Table 4 

HETG Spectral Parameters of 2 Temperature APED fits 


Name 

Nh 

Ti 

T2 

EMi 

EM2 

Fx 

Lx 


(1) 

[MK] 

[MK] 

(2) 

(2) 

(3) 

(4) 

MTOri 

i.JJ Q J4 

11 00 

11.UU 0 47 

39.24 };| 

1 2 

0.1 

6.6 ®;2 

20.2 {j 

8.4 

LQOri 

0.00 0;^® 

Q on 0.27 
0.27 

39.10 li 

0.1 

3.4il:j 

13.8®:^ 

5.0 

Par 1842 

0.51 O;?? 

1 1 ng 0.54 

37.50 1^2 

0 7 

0.1 

46.0 

3.1 

V348 Ori 

0.83 °;j« 

1 1 A4 0.46 

45.2311 

0.9®:® 

3 3 

0.1 

35.4 |j 

4.6 

VI229 Ori 
VI399 Ori 

0.84|j;20 
0.72 [!:j« 

10 34 0.52 
lU .:)^ 0.88 

0.37 

36.41 \ll 

37.96 

1-0 

3 0 

3 0 

0.1 

25.5 };® 
31.8 \l 

3.9 

4.3 


(1) 10^' cm ^ (2) 10^"^ cm ^ (3) 10® erg cm ^ s ' (4) 10^1 erg s ' 


Table 5 

Line Eluxes of Major Lines used for Diagnostics 


Line 

MTOri 

LQOri 

Par 1842 

V348 Ori 

V229 Ori 

VI399 Ori 


(1) 

(1) 

(1) 

(1) 

(1) 

(1) 

Ar XVIII La 

0.80 

0.34 0:i 

0 19 

0.12 

0.04®:!® 

0.03®:®! 

0 19 ®12 

Ar XVII r 

1.08 0:11 

0.28 0:10 

0.22 

0 22 

0.21 ii:l? 

0.31 

S XVI La 

0.75 U:!® 

0.43 0:1 

A 07 0.30 

0.43 ®:i 

1.48 0:23 

0.38 0:i 

SXVr 

1 17 

0.43 

0.58 0:21 

A QQ 0.38 
0.38 

0.48 ®:11 

0.38 ®:22 

U.J/ 0.29 

SXVf 

A 7A 0.35 

U. /U 0 40 

A ^4 0.50 
0.04 

A 43 0.26 
0.26 

0.30 °:i 

0.58 0:1® 

A 29 0.22 

Si XIV La 

0.29 

1 70 0.22 
0.23 

A AO 0.21 
0.21 

1 26 ^-25 

1 14 0.21 

0.22 

1 43 0.22 

Si XIII r 

1 43 0.28 

^ 0.29 

U.DJ 0.21 

A 0.20 

0.20 

0.72 0:28 

1 A4 0.23 
0.25 

1 A2 0.20 

0.22 

Si XIII f 

0.25 

0.52 :1? 

0 32 

0J7 

0.43 p2 

U.Dl 0.23 

0 42 01® 

U.H.^ 0.17 


Mg XII La 

1 07 0.28 
0.30 

0.80 0:1? 

0 47 ®i* 
^2.^' 0.18 

A 92 0.20 
0.22 

1 10 

Mg XI r 

0 27 

0.24 

0 32 

U.ZD 0.17 

A 77 0.25 
' 0.28 

U.JJ 0.23 

0 71 
^ 0.22 

Mg XI i 

0 76 

UOJ 0.32 

A AO 0.00 
0.00 

U.i:) 0.24 

0 26 
u.zD 0.21 

0 42 

Mg XI f 

0.20 0:20 

0 00 00® 

0.33 ®:®® 

A 29 0.20 

0.24 

0.27 0:1? 

0.47®:!! 

NeXLa 

24.121:1 

20.06 fll 

8.86 j:°l 

1 1 39 4.48 
1.38 

12.64 I:!' 

8.331:11 

NeIXr 

4.33 2-}® 

4.201:00 

2.64 ®:«j 

7 77 0.98 
1.15 

2 40 

0.95 

117 0.60 

0.73 

NeIXi 

0.83 0:®2 

0.84 0:12 

0.85 ®:®1 

0.92 

1 A7 0.56 
i.U/ 0.71 

0.54 

NeIXf 

A 40 0.83 
0.97 

1 77 0.74 
0.88 

1 74 0.72 
0.72 

1 1 1 0.68 
0.88 

2.07 in 

1 AO 0.55 

0.75 

O VIII La 

16.281:11 

29.59 

5 21 

2.46 

Q.QO 

2.381:11 

0 70 2.87 

3 56 

4 94 T79 

6.07 

O Vllr 

6.011®2 

0.20 O:?® 

7 65 

4.73 

2 55 

6.92 

0.00 ®:i 

O Vlli 

3 01 

0.201:®® 

4.141:1® 

6 35 

8.69 

12.81 lo'?o 

2 07 2.17 

5.30 

O Vllf 

7 1 7 ®-27 

3.70 I'll 

5.01 III 

3.54 l:?l 

0 20 

13.19 ?o'^6 

N VII La 

6 42 

7.01 

4 QQ 3.63 
5.82 

0 00 ®®® 

0.00 

0.00®:®® 

6 17 ®'22 

' 8.05 

8.36 ®:21 


(1) 10 ® photons cm ^ s ' 
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dents a, „ ai, and a?, were determin ed from APED and are 
listed in Huenemoerder et’aP ( 1201 3h assuming abundances of 
lAnders & Grevessel(ll989l) . 

The temperature insensitive abundance ratios we calculate 
from this procedure are listed in Table |6] and plotted in Eig. 0 
There are more recent and update d distri butions such as the on e 
from iGrevesse & Sauval d 19981) and lAsplund et alj (|2009|) . 
The differences, however, are well within the uncertainties of 
the X-ray abundance ratios an d we did not re-scale the values 
to the most recent distribution (lAsplund et al.ll206^ . 



Fig. 5.— Measured abundance ratios for Orion sources. The symbols mean: 
circle(Ne/Mg); cross (Mg/Si); square (Si/S); tr angle (S/0); diamond (Ne /0); 
star (S/Ar). Added are the values from the sun JGrevesse & Sauvalll 19981) and 
the ISM jWilms, Allen & McCradl2000l) 

The average value of these ratios for each ratio type was 
the basis to construct a set of abundances and associated un¬ 
certainties. These were used to constrain the abundance pa¬ 
rameters for O, Ne, Mg, Si, S, and Ar during the APED 
plasma fits below. We also compare these ratios to the val¬ 
ues e xpected from the Sun dOrevesse & Sauvall[l998h and the 
ISM dWilms. Allen & McCravll2000l) . We observe that some 
values such as Mg/Si or Ne/Mg or Si/S align fairly well with 
the Sun and the ISM, S/O or Ne/O or S/Ar are significantly dif¬ 
ferent indicating anomlies with the oxygen and/or neon abun¬ 
dances. In our case O Vllwas only poorly detected, however. 
The abundance of Ee was not included in this procedure for the 
lack of K-shell lines and ratio coefficients. Here we only can 
adjust the abundance during the APED fits. 

4.2. R-Ratios from He-like Triplets 

In collisionally ionized plasmas the He-like line triplets are 
sensitive to density and external UV fields. He-like line triplets 
consist of resonance (r), intercombination (i), and forbidden (f) 
line components. At increasing collisional plasma densities the 
meta-stable f-line component depopulates into the i-line com¬ 
ponents. A similar effect appears by the influence of far-UV 
photons matching the energy difference between the f- and i- 
lines. 

In the absence of a strong UV radiation field, the flux 
ratio {R = f /i) of the two lines is sensitive to the density 
of the plasma. Coronal plasmas are considered low den¬ 
sity plasmas with electron/ion densities below the order of 
10'* c m~^, while plasmas in accretion columns can reach 10*"* 
cm“^ (iKastner et al.ll200^ . He-like triplets of O, Ne, and Mg 
ions are sensitive to this range. In general, R-ratios signifi¬ 


cantly below 4.5 indicate densities significantly exceeding 10*" 
cm“^ (iTesta, Drake & PeresI l2004l) . Eor neon, ratios signifi¬ 
cantly below 3.8 indicate densities significantly exceeding 10** 
cm“". And last, but not least, for magnesium, ratios signifi¬ 
cantly below 2.8 indicate densities significantly exceeding 10*^ 
cm“". However, even though many R-ratios in our sources are 
below these limits, in nearly all but a few cases the result is not 
very significant because of the low statistic in the lines. At this 
stage the uncertainties do not allow us to discern from the low 
density limit. 





Wavelength (A) Wavelength (A) 


Fig. 6.— APED plasma fits of the Ne IXtiplet regions in the six Orion 
sources . 

Figure |6] shows the Ne IX triplet region for all sources. Of all 
the lower Z triplets the Ne IX triplets appear more prominently 
in our spectra than the Mg XI and O VII triplets, also likely due 
to a neon overabundance. Weak He-like line strengths are likely 
due to the high plasma temperatures these spectra, which favor 
resonance line strengths. This can be seen in Figure |6l where 
in all sources the Ne IX resonance line dominates the bandpass. 
The Ne IX bandpass in Figure |6] was fit with the lower tem¬ 
perature APED plasma component and the results shows that a 
plasma of low density (< 10** cm“^) fits the region well within 
the limited statistics. 

4.3. Column Densities 

The broadband plasma fits produce unusually low column 
densities, some with values well below expectations from the 
COUP analysis, suggested values from observed optical ex¬ 
tinction, and sugg estions from previous analysis of the m assive 
stars in the ONC (ISchulz et akOOPlIlGagne et al.ll2005h which 
suggest columns of at least 2.9 x 10^ * cm~^. Columns can be 
estimated from optical extinction via (iPredehl & Schmit3l 19951) 
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Table 6 

Temperature insensitive line ratio with uncertainties 


Name 

Ne/Mg 

Mg/Si 

Si/S 

S/0 

Ne/0 

S/Ar 

MT Ori 

10.2±4.8 

0.7±0.2 

0.6±0.2 

0.4±1.3 

11.3±7.5 

0.5±0.2 

LQ Ori 

12.7±6.4 

0.6±0.2 

0.8±0.4 

3.0±2.7 

2.0±1.8 

0.8±0.5 

Par 1842 

7.4±3.3 

0.4±0.1 

0.4±0.2 

- 

- 

0.9±0.6 

V348 Ori 

4.5±1.6 

0.4±0.3 

0.9±0.5 

- 

- 

1.4±1.0 

VI229 Ori 

5.6±1.8 

0.5±0.2 

0.5±0.3 

O.liO.l 

14.9±13.4 

2.4±2.2 

VI399 Ori 

3.0±1.1 

0.6±0.2 

1.0±0.4 

- 

- 

0.9±0.6 

4 ' Ori E 

1.7±0.3 

l.OiO.l 

1.2±0.2 

1.4±1.0 

4.9±2.8 


Sun 

3.16 

1.10 

1.91 

0.02 

0.14 

4.12 

ISM 

3.46 

1.35 

1.51 

0.03 

0.18 

4.79 



from [10^' cm-^] 

Fig. 7.— The column densities as predicted from optical extinction Ay ver¬ 
sus the final columns found in the plasma model fits. The bold line indicates 
the expectation of the two values being the same. 


Av = 0,56 xN„-0.23, 


( 2 ) 


where Ay is the extinction (see Table|2]l and Nh is the X-ray 
column in units of 10^^ cm“^. Figure|2]shows the columns from 
optical extinction Ay versus the measured X-ray columns. Here 
we used the columns from the APED model fits. These columns 
are significantly lower than the one expected from extinction. 

We have also determined line-of-sight absorption from line 
fluxes of approximately known intrinsic ratios. Once emission 
lines are from a single ion, the ratios are determined by the 
atomic data, and to a weaker extent, on the plasma temperature 
distribution. If these intrinsic ratios are known, then the ob¬ 
served ratios are determined by the absorption function, which 
then can be inverted to give the column density. We used the H- 
like Lyman-like series with O VIII a (A18.973) and /3 (A16.006) 
lines. The Ly-^ like line, however, is blended with Fe XVIII 
A16.004 (lower and upper transition configurations of - 

2p^). 

Using another Fe XVIII line at 14. 208 (a doublet; 2p‘^3d - 
2p^), we can break the degeneracy (iTesta et akl |2004 l2007h . 
We have three measured fluxes comprised of 4 lines, with two 
known ratios (an Fe XVIII ratio, and the O VIII ratio). The 
other unknown is the absorbing column. This constitutes five 
equations with five unknowns, so we can solve for Ah. The 
abundances are not a parameter since we are ultimately using 
the ratios of identical ions. 

The system is non-linear (since the absorption function is 
exponential), so we solved this numerically. We used the 
“phabs” function f rom XSPEC l Arnaudll 19961) and the emis- 
sivities of AtomDB (iFoster et al.ll2Q12 1. We show an example 
fit to the relevant MT Ori spectral region in Figure [8] 



Wavelength (A) 


Fig. 8.— The MT Ori region from which we derived the absorbing column 
from emission line ratios. The fit function was a constant plus four Gaussian 
functions. The lines of interest are Fe XVIII A14.208, O VIII A18.967, and 
the unresolved blend of Fe XVIII A16.004 and 0 VIIIAI 6 .OO 6 . We also in¬ 
cluded a resolved blend of Fe XVIII A16.071 in order to get a better flux of the 
unresolved blend. 


Given the weakness of the 14.208 A line, we performed a 
Monte-Carlo analysis in which we perturbed the line fluxes 
according to their uncertainties (assuming Gaussian distribu¬ 
tions). Distributions of the column density were then obtained 
for 10"^ realizations. We did this for the PMS stars in this sample 
as well as for the more massive ONC star s, and for two bright, 
unabs orbed objects, crGem and HR 1099 (iHuenemoerder et^ 
l2013h . We also added a distributon for the sum of the Orion 
PMS sample stars. The Orion PMS stars consistently yield a 
most probable value of Ah 0.2 x 10^^cm“^. The “control” 
stars (crGem and HR 1099) yielded much lower values of about 
0.03-0.05 which is consistent as they are expected to be con¬ 
sistent with almost 0.0. Figure |9] shows the distributions ob¬ 
tained for the set of Orion stars analyzed. The peaks of line 
ratio distribution is mostly consistent with the columns from 
the broadband analysis, LQ Ori appears a little bit higher but 
the 90% uncertainties align well with the propabilites just be¬ 
low the peak. We conclude that the columns in the Orion stars 
in this sample are indeed very low. 
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0.01 0.02 0.05 0.1 0.2 0.5 

N„ [10^2 cm-=] 


Fig. 9.— The fits to the collection of stars, as determined from 10^ Monte- 
Carlo realizations for each star. 

4.4. APED Fits: Temperatures and Abundance 

Broadband fits of all spectra were performed in several steps. 
The first step involved a restricted fit using a first set of abun¬ 
dance settings obtained from the temperature insensitive line 
ratios. The uncertainties in the set are large and allowed the 
fit considerable freedom. In consecutive steps we manually 
adjusted single elements during the fits in order to minimize 
residuals which appeared in the He-like and H-like lines. Eor 
the final distibution we calculated 90% uncertainty limits for all 
values. The final abundance distribution fo r the fits then yields 
the fo llowing values with respect to solar (I Anders & Grevessa 
119891) : O (0.195±0.38), Ne (0.718±0.23), Mg (0.085±0.034), 
Si (0.153±0.07), S (0.225±0.120), Ar (0.524±0.220), Ee 
(0.058±0.012). 

Multi temperature components did not produce statistically 
superior solutions with respect to two temperature models. The 
final two temperature fits are shown in Fig.[3and Fig.|4]and the 
parameters are shown in Tab. |4] The scatter on most spectra 
is very large above 20 A and for the final APED fits we limited 
the bandpass to 1.8 and 20 A. The data bins were grouped to the 
size of 0.04 A in the brightest and 0.08 A in the weakest source 
for display purposes. For the fits, each spectrum had 941 inde¬ 
pendent data bins, each bin had the size of one MEG resolution 
element (0.021 A). 

The resulting X-ray plasma properties we measure in our 
sample of six sources, i.e. emissivities and temperatures, gen- 
erally agree well wi th the overall COUP sample trends in 
iPreibisch et al.l ( l2005l) . This is, h owever, not true for the source 
parameters measured in COUP (iGetman et al.ll200'5l) . A com¬ 
parison of Table [3 with Table |4] shows that fluxes and emissiv¬ 
ities are overestimated and temperatures significantly lower in 
the COUP data, which likely adds to the scatter in the COUP 
results. 

5. RESULTS AND DISCUSSION 
5.1. Flux and Fuminosity 

Even though the limited number of sources in this sample 
does not allow for any statistical trend analysis, we can point 
out some differences and similarities between these sources. 
There are some differences in the overall luminosity, the weak¬ 
est source Par 1842 has 3.1x10^' erg s“* , the brightest source 
MT Ori has 8.4x 10^' erg s“* . The spread increases when we 
compare the surface flux, i.e. the X-ray luminosity divided by 
the stellar surface area. Figure [Tol shows the surfac e flux with 
respect to the listed ages from iDa Rio et alJ (1201 Ol) . which in¬ 


dicates that in the age range up to 10 MK we observe a sub¬ 
stantial increase in coronal activity. This trend is implicitely 
reflected when estimating coronal loop sizes, in this case using 
em ^ p^r^i, where em is the volume emission measure, p a coro¬ 
nal density of 10^' cm“^, and r; the loop size. This is a very sim¬ 
plistic estimate, we here also just assume a Ailing factor of unity 
and similar densities, but it shows that these sizes also increase 
with age. The fact that surface normalized X-ray fluxes increase 
with age for CTTS in the ONC has been repo rted by COUP, but 
for lower mass stars dPreibisch & Feigelsonll200^ up to 1 Mq 
but with no significant effect for 1-2 Mq stars. Here stars 
range from 0.7 to 2.3 Mq and the effect is highly significant. 
One should take the case of LQ Ori with caution, however, its 
age determination of < 0.1 Myr is unreal istically low and no t 
supported by models, but based on the data lDa Rio et al.l (1201 Oh 
it is still likely younger than MT Ori. 

5.2. Column Densities and Extinction 

Except for LQ Ori all sources in the sample are affected by 
substantial amounts of extinction (Tab. |2]l, which is not un- 
us ual and eyen expec t ed tow ards the ONC. In ROSAT stud¬ 
ies [ftedehT^^chmit^ (|l995t) inyestigated the relation of X-ray 
column density and optical extinction and found a quite robust 
relationship (Eq. |2]) specifically yalid for columns aboye 10^' 
cm“^. Howeyer it was deyeloped for long range line of sights 
with gas and dust distributed oyer yast distances. In the case of 
Orion we haye the situation that extinction is local to the dust 
in the ONC and its sources, while the X-ray column reflects ab¬ 
sorption towards Orion and is likely gas dominated. In other 
words, there must be a different offset than in Eq.|2]due to local 
dust, which is not obseryed in the X-ray column. Such an offset 
in Ay is clearly yisible in Fig. 0 

Howeyer, such interpretations haye to be yiewed with cau¬ 
tion. In star forming regions the extinction is indeed local and 
the extinction law might differ from the properties in the in¬ 
terstellar medium (ISM). In particular, the properties of dust 
grains, which cause the optical extinction, are expected to 
change as the star forming regions eyolyes. In the initial cool 
starless cores temperatures are much lower and densities much 
higher than in the ISM, fayoring dust formation, while later 
the intense radiation from the newly-formed stars will eyapo- 
rate smaller dust grains. Obseryationally, one way to probe this 
regime is to compare the X-ray absorption Nh with the optical 
extinction Ay. X-ray absorption is caused by the total column 
density (gas and dust) of heayy ions and is expressed as the 
equiyalent hydrogen column density, assuming some standard 
set of abundances. On the other hand, the optical extinction is 
dominated by the dust colu mn density and also depends on the 
dust grain size distribution. IVuong et akl (l2003h compared Nh 
and Ay for six star forming regions, including the ONC. For the 
p Oph star forming region they And an Nh /Ay significantly be¬ 
low the ISM yalue and they interpret this as a sign that the cloud 
material in p Oph has a lower metal abundance than the ISM, 
consistent with recent solar abundance measurements. Alterna- 
tiyely, grain growth can increase the amount of e xtinction per 
unit m a ss until the grains re ach about 1 pm in size (lOrmel et al.l 
1201 ih . IVuong et akl (l2003h haye only a few data points with 
a large scatter in the ONC, which are roughly consistent with 
an ISM-like Nn/Ay ratio, but in our data (Fig. 5) we see that 
all targets haye an N n/Ay ratio significa ntly below the ISM 
yalue, in line with the lVuong et akl (l2003h result for p Oph. In 
fact, low Nh/A y yalues are seen in seyeral low-mass star form- 
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ing regions, e.g. NGC 1333 ( Winston et al.ll2010h . and IRAC 
20050+2720 dGiintheret al.ll2012h . Confusingly, pr evious stud¬ 
ies of high-mass st ar forming regions (RCW 38: IWolk et ^ 
(l2006l) : RCW 108: IWolk et alJ (l2008h ') tend to hnd ISM-like 


Nh/A y ratios. It should be noted, however, that the studies cited 
above all use CCD spectra to measure Nh and that there often 
is an ambiguity between the amount of cool plasma and the ab¬ 
sorbing column density, such that a high column density and a 
large amount of cool plasma or a low column density with no 
absorption provide and equally good ht. In our study we break 
this dege neracy through the use of high-resolution grating spec¬ 
troscopy. IWolk et ^ (l2008h used a similar approach for older, 
but still accreting CTTS in low-mass star forming regions and 
found consistently high Nh/A y ratios, indicating gas-rich ma¬ 
terial. They speculated that some of their sources could be seen 
through the accretion column, which would be dust-depleted 
due to the stellar irradiation. 

We are left with a complex picture where the dust content, 
grain growth, and abundance in the circumstellar matter and 
the ISM all influence the observed Nn/Ay ratio and their re¬ 
spective influence in the sightline to a specific target cannot be 
disentangled. For our six targets in the ONC, the total extinc¬ 
tion is low, so there cannot be much cloud material in the line 
of sight. Given that the objects are all very young, some cir¬ 
cumstellar matter with evolved dust grains and possibly a low 
metal abundance can explain the observed A///Av ratios. 
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Fig. 10.— The stellar surface X-ray flux versus stellar age. 


5.3. Plasma Densities and Accretion 

The fact that we cannot use observed R-ratios because of 
our limiting statistics at this stage is unfortunate, but one also 
has to put their value for our star sample into perspective. At 
these high plasma temperatures even for the lower tempera¬ 
ture component, the He-like triplets are not only difficult to ob¬ 
serve, but we would not see much contribution from accretion 
shock plasmas in these sources. These stars have simply not 
contracted enough to support sufficiently high infall velocities, 
which range from as low as 150 km s“' in LQ Ori to about 350 
km s“^ in Par 1842. For our pred ictions we use the v a lues i n 
Tableland a simple model from ICalvet & Gullbringl (Il998h . 
For example, in order to ionize the plasma up to Mg XI temper¬ 
atures above about 2 MK are necessary, which requires infall 
velocities above 800 km s“'. This cannot be the case here and 
we would not expect signatures from accretion activity. 

Even though our hts of the Ne IX regions are consis¬ 
tent with a low density plasma, we are statistically lim¬ 
ited to rule out higher density plasmas. These might be 


observed in specihc coronal environments at the basis of 
strong coronal lo ops involving small emitting volumes and 
high temperatures d Sanz-Fo rcada, Brickhouse & Dupr^l2003t 
iTesta, Drake & Pere*^ l2004fi or coronal heating of dense 
and fla red inner a cc retion disk regions (Drake & Orlandol 
|2010t Drake et al. _ 20091 iDullemond & Monniej 20101 : 


iDullemond et al. 


_ 2OO7I) . Magnetic conhnement as outlined 

by Gagne et al.l ( 2005T might be a possibility in these type of 
interactions as well. 

In the following subsections we dicuss how these issues fare 
with individual sources and how these properties compare with 
other PMS stars and results from the COUP sample. 


5.4. Individual Sources 

MT Ori is the X-ray brightest low-mass PMS star in our 
sample and with a surf ace temperature of 4 600 K currently a 
K2 - K4 spectral type (iDa Rio et alJl2010h . The latter study 
determines a contribution of about 13% of the total luminosity 
from accretion. At an age of about 0.5 Myr (see Tab. |2]), one 
might expect more from accretion in a 2 Mq CTTS. The ratio 
of the emission measures of the two temperature components 
yields a fraction of 18% for the lower temperature component. 
An accretion fraction would be consistent with the low optical 
accretion luminosity, however, the estimated infall velocity of 
250 km s“' is too low to allow for K-shell ionization of matter 
with atomic number larger than O VII. We should thus not ex¬ 
pect low R-ratios from an infall accretion shock for Ne IX and 
Mg XI. 

The total X-ray luminosity is very high leading to a log 
Ljc/Lhoi of abo ut -2.7 , abov e the saturation limit of -3 for fast 
rotators (e.g. IVilhul (11984 : IVilhu & Waited (Il987h ') and well 
above the median value of -3.7 for TTS with unknown periods 
in COUP data. However, the surface X-ray flux is the 2nd low¬ 
est in the sample, due to the fact that the stellar surface area 
is the largest in the sample. The lightcurve shows the source 
persistently bright at levels varying by a factor of 8. From the 
average total emissivity we project average coronal loop sizes 
of about 10% of the stellar radius assuming coronal densities 
of 10" cm“^. iGetman et alJ (l2008h found a typical flare in the 
COUP dataset and estimate loop sizes of about 1.1x10*' cm, 
which is only 30% larger. Note, the estimated loop sizes use a 
model by [Reale & Micelal(ll998h . which depends on the square 
root of the observed flare temperature. However, because of 
pileup problems (see Sect. I4.41 i this temperature is about a fac¬ 
tor two lower in the COUP dataset and the true sizes are larger 
by a factor s/2.3. 

LQ Ori is the second X-ray brightest low-m ass PMS star 
in our sample and classified as a K2 V dwarf. iDa Rio et al.l 
(1201 Ol) And less than 1% of the optical luminosity due to ac¬ 
cretion, which for a 0.6 Mq star at an age of the order of 0.1 
Myrs is highly unusual. Ages this low are not supported by 
most model calculations and we expect LQ Ori to be somewhat 
older. The fraction of the low temperature component is almost 
40% indicating the potential for some accretion to produce X- 
rays. However, the case for X-rays from accretion in LQ Ori 
is even worse than for MT Ori as predicted infall velocities do 
not even favor much ionization up to O VII. LQ Ori also does 
not appear to be a particularly dusty environment, extinction 
is extremely low in this star and this is backed up by a very 
low column density allowing to record even a strong N VII Lq 
line, log LxILhoi is slightly lower than the one in MT Ori, but 
still well above -3. But similar to MT Ori we estimate an av- 
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erage coronal loop size of about 10% of the stellar radius. A 
inco mplete flare in the CO UP data provided about 2.5 times the 
size (iGetman et al.ll2008h . 

Par 1842 is the oldest star in the sample according to 
iDa Rio et alJ (1201 Ol) with an age of 10 Myr. Its determined lu¬ 
minosity fraction from accretion is the highest in the sample. 
The low temperature component produces about 30% of the to¬ 
tal emission measure, and while models predict infall velocities 
of the order of 350 km s“' allowing for both, ionizations up to 
O VII and Ne IX but not to Mg XI are feasible. The X-ray lumi¬ 
nosity is the lowest of all stars, but its surface flux the highest, 
indicating that this star is actually coronally most active. Con¬ 
sequently from the total emission measure we estimate average 
loop sizes of over 20% of the stellar radius. One flare in the 
COUP data, which was classified as a slow rise event, gave a 
loop limit over 10 time s larger (5x 10*^ cm ) indicating that this 
flare was a giant event (iGetman et alj|200^ . 

V348 Ori, V1229 Ori, V1399 Ori have very similar op¬ 
tical pro perties in terms of surface temperature, mass, radius, 
and age (iDa Rio et alj|2010|) . Overall this is also true for their 
X-ray properties. The optical luminosity from accretion (see 
Table |2]i, however, varies by about a factor two between VI229 
Ori, V348 Ori and VI399 Ori likely indicating that accretion 
rates vary by factors independent of system parameters. R- 
ratios in all three sources show reduced values and the presence 
of cooler dense plasma from accretion cannot be ruled out. The 
case here is similar to the one for Par 1842, ionizations in the ac¬ 
cretion shock plasma might favor O VII or Ne IX but not Mg XI. 
log is also above -3 and average coronal loop sizes are 

of the order of 20% o f the stellar radius. N o flares were detected 
in the COUP sample dGetman et al.ll200^ . 


5.5. Comparison with other T Tauri stars 

There are few detailed studies of young CTTS exhibit¬ 
ing dominant coronal properties other than the Orion stud¬ 
ies. High resolution X-ray studies dealing with young PMS 
stars at various evo lutionary stages have been conducted by 
lAudard et all dlool for V987 Tau, RY Tau, and LkCa 21, or 
by [Robrade & Schmid (l2006h for the four CTTS BP Tau, CR 
Cha, SU Aur, and TW Hya. Perhaps the most prominent non- 
Orion example from these samples is SU Aur. It was found 
that coronal emission is the dominant source of X-ray activ¬ 
ity in this star with most emissions between 20 - 50 MK and 
the coolest component at aroun d 8 MK . Due t o absorption, the 
O VII triplet is not observed. IShuklal (l2009t) And R-ratios of 
Mg and Ne that limit plasma densities to less than 10'^ cm“^ 
and conclude a dominant coronal nature. This is very similar to 
what we find for the ONC sources, except for the fact that ONC 
sources appear significantly more luminous. Consequently SU 
Aur appears below the Orion source data points in Fig. [TO] 

Most PMS coronal sources are classified as weak lined T 
Tauri Stars (WTTS). Here the IR-excess from a prominent ac¬ 
cretion disk has vanished and we do not expect any accre¬ 
tion activity. X-ray emissions are considered to entirely be of 
some form of coronal origi n. Prime examples are HD 98800 
in the TW Hya association (iKastner et al.ll2004li a nd V987 Tau 
in the Taurus-Auriga region ( Strassmeier & Ric3ll998h . These 
non-flaring states of WTTS are well described by either a one 
temperature component or a weak hot and dominant very hot 
component. The case of V987 Tau is of interest because it s 
age is 2 -3 Myr (Furlan et ^120061 : [Pamiani & Micelalll995h . 


[Andrews & WilliamsI (l2005h found a disk mass of < 0.0004 


Mfn in dicating that accretion has ceased in this system. IShuklal 
(l2009l) estimate R-ratios of Mg and Ne that limit plasma densi¬ 
ties to less than lO'^ cm“^, however for a very weak cool com¬ 
ponent. In Fig. [To] V9 87 Tau appears even further below th e 
Orion stars and SU Aur (ISiwak et alj201 ll: lAudard et al.l200^ . 
which is likely a consequence of its much earlier spectral type. 

5.6. Abundances in the ONC 

The ONC as an ensemble of coeval stars allows to study de¬ 
tails in the chemical composition of these stars and the clus¬ 
ter. Coronal chemical compositions are often diffe rent from 
stella r photospher ic and interstellar compositions dFeldmanl 
1 19921: lTbstall201^ . Most of these differences relate to deple¬ 
tion effects with respect to the first ionization potential (FIP). 
The solar corona has more elements with low FIPs such as 
Mg, Si, Fe, while coronally active stars seem to have more 
higher FIP elements such as O, Ne, Ar. Figure [TT] com¬ 
pares the average abundance ratios of th e Orion stars with 
various studies. Solar pho tospheric ratios (iGrevesse & Sauv^ 
ll998L[Asplund et ani2009l) seem to be quite in agreement with 
abunda nces from B-stars (Asplund et al.ll200^ and B-stars in 
Orion (iNieva & Simon-Dl^ 2011 ). A study of the upper so¬ 
lar atmosphere at 10® K seems to also differ significantly from 
the Orion ratios. The ratios for the Orion stars seem to agree 
with most of these studies only for the Ne/Mg, Mg/Si, and 
Si/S ratios, while significant differences can be found with re¬ 
spect to the S/O, Ne/O, and S/Ar ratios. This indicates that 
there are likely discrepancies with at least the neon and oxy¬ 
gen abundance. Enhanced neon has bee n observed in X-ray 
spectra of other stars than the Sun (see iDrake & Testal l2005l 
and reference therein). Besides enhanced neon, the very high 
Ne/O ratio could also indicate depletion of oxygen. Given that 
column densities are not extraordinarily high and that we ob¬ 
serve the N VII line in some cases shows that there could be a 
depletion of oxgyen in the Orion spectra. This can have sev¬ 
eral reasons ranging from the fact that spectra at these high 
temperatures mostly generate H-like O VIII lines, but also 
due to the large extinction values, which indictate that these 
CTTS are likely quite dusty and oxyge n is depleted into dust 
grains (iDrake, Testa & Hartmannir2005l) . 
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Figure[T2]shows a similar comparison but now involving stel¬ 
lar active coronae. He re we used two stu dies, one for a set of 
fast rotating M dwarfs (iLiefke et al.||2008l) and two active coro¬ 
nal K stars (iHuenemoerder et al.ll2013l) 7Even though for some 
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ratios the uncertainties are still large and there are large varia¬ 
tions within Orion sources there is a very good agreement with 
ratios of coronally active sources. Even the Ne/O ratio with 
its quite large uncertainties seems more consistent with coronal 
souce s dHuenemoerder et al.l 120131: Drake. Testa & HartmannI 
l2005h than with solar values (lAsnlund et al.ll2009ll . The com- 
parison of Fig. [TT] and Fig [12] thus highlights that elemental 
abundances in stellar photospheres differ from stellar coronae. 
This difference in abundances is still an unsolved problem. 

The measured abundance ratios also show that they are in¬ 
sensitive to a range of very high coronal temperatures but 
sensitive to a more global temperature range. However, ex¬ 
actly this difference between photospheric and coronal abun¬ 
dances can work i n favor as a d i agnostic probe for coro¬ 
nal structures (see Lam ing 1201^ 

Sanz-Forcada. Brickhouse & Dupree 


Laming & Hwang! l2009t 
20031 testa2015') In our 


Orion sample we are still far from doing so given the current 
statistical limitations in the X-ray line fluxes. 

Abundances of iron could not be included in the line ratio 
study but from the spectral plasma fits they appear consistently 
low in all Orion source s. Low Fe abundances are also found 
in active stellar coronae dHuenemoerder et '^120 1 3h but also in 
X-ray spectra of the Orio n Trapesium such as the in termedi- 
ate mass binary 0' Ori E dHuenemoerderet^L 2013h and th e 
massive stars 0' Ori C and 6^ Ori A ( Schulz et alJl2003Ll2006h . 
However, in all these cases that observed plasma is likely of also 
coron al origin. Other studies of Orion dNieva & Simon-DfazI 
1201 Ih and the solar neighborhood do not support such an iron 
abundance deficit, pointing once again to chemical fractiona¬ 
tion effects in the outer stellar atmosphere. 
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Fig. 12.— The average abundance ratios for our Orion sources (circles) i n 
comparison with various M dwarfs with active coronae ILiefke et alj|200^ . 
The symbols mean: cross (YY Gem); square (AU Mic); trangle (EV Lac); di¬ 
amond (AD Leo); star (Proxima Cen); filled circle (a Gem); filled square (HR 
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6. CONCLUSIONS 

We observed some of the youngest CTTS in the ONC with 
ages between 0.1 and 10 Myr. All stars are very bright in X- 
rays with luminosities significantly above 10^' erg s“* and ra¬ 
tios of X-ray to bolometric luminosity of the order of -2.8. The 
X-ray surface flux increases significantly with young age even 
for stellar masses above 1.5 Mq together with average coronal 
loop sizes indicating an increase in coronal activity during the 
CTTS phase. The column densities towards the X-ray sources 
are lower than expected from measured extinction values and 
lower than previously determined for the central massive star 
0* Ori C, likely due to high local gas-to-dust ratios. The X-ray 
spectra are too hot to allow for significant R-ratio determina¬ 
tions with current statistics, but on average the He-like triplets 
are consistent with typical densities of coronal plasmas. Aver¬ 
age abundance ratios are also consistent with coronally active 
stars underlining the coronal nature of the observed X-ray emis¬ 
sion. The preliminary results of young coronally active PMS 
stars in this analysis demonstrate the incredible potential of ob¬ 
serving the ONC stars with the Chandra HETG to study details 
of the coronal evolution of CTTS. 
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